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ABSTRACT 

Recent high resolution spectroscopic analysis of nearby FGK stars suggests 
that a high C/O ratio of greater than 0.8, or even 1.0, is relatively common. Two 
published catalogs find C/0> 0.8 in 25-30% of systems, and C/0> 1.0 in ~ 6- 
10%. It has been suggested that in protoplanetary disks with C/0> 0.8 that the 
condensation pathways to refractory solids will differ from what occurred in our 
solar system, where C/0= 0.55. The carbon-rich disks are calculated to make 
carbon-dominated rocky planets, rather than oxygen-dominated ones. Here we 
suggest that the derived stellar C/O ratios are overestimated. One constraint 
on the frequency of high C/O is the relative paucity of carbon dwarfs stars 
(10~ 3 — 10~ 5 ) found in large samples of low mass stars. We suggest reasons for 
this overestimation, including a high C/O ratio for the solar atmosphere model 
used for differential abundance analysis, the treatment of a Ni blend that affects 
the O abundance, and limitations of one-dimensional LTE stellar atmosphere 
models. Furthermore, from the estimated errors on the measured stellar C/O 
ratios, we find that the significance of the high C/O tail is weakened, with a 
true measured fraction of C/0> 0.8 in 10-15% of stars, and C/0> 1.0 in 1-5%, 
although these are still likely overestimates. We suggest that infrared T-dwarf 
spectra could show how common high C/O is in the stellar neighborhood, as the 
chemistry and spectra of such objects would differ compared to those with solar- 
like abundances. While possible at C/0> 0.8, we expect that carbon-dominated 
rocky planets are rarer than others have suggested. 
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Introduction 



1.1. The Composition of Stars and Planets 

The determination of the abundances of atoms in the atmospheres of stars is an es- 
sential element of modern astronomy. Recently, tremendous work has occurred on under- 
standing the relationship between planets and the a bundances of pla net-hosting and non- 
planet-hosting stars. Since the pioneering work of iGonzalezI (119971 ). many investigators 
have worked to understand connections between stellar abundances and the observed fre- 



uencv (pantos et a 



Guillot et al. 



2006 



2004; 



Fischer fc Valenti 



2005: 



Johnson et al.l 120101 ) and composition 



Burrows et al.ll2007t iMiller fc Fortnevlboilh of planets. 



Our solar system is one realization of the complex planet formation process. The raw 
materials that made up the Sun and solar nebula, through a process of condensation, grain 
growth, and accumulation, gave rise to four rocky planets in our inner solar system that are 
predominantly composed of Mg-Si-O-bearing rocks and Fe-Ni metals. In other solar systems, 
with parent star disks with other abundances, a different selection of refractory materials, or 
in different relative abundances, surely occur. For instance, if a nebula's carbon-to-oxygen 
(C/O) ratio is >0.8, condensation pathways can change dramatically , leadi ng to carbon- 
dominated rocky planets, as recently discussed in detail by iBond et al.l ( 12010I ). 



There h ad been prior in termittent inte rested in carbon-dominated pla nets in the past 
decade, from Gaidos ( 2000 ) . lLodders ( 2004 ) . and Kuchner fc Seager ( 2005 ). to name three 
examples. In particular, iGaidod ( 120001 ) discussed different formation scenarios for giant 
planet cores and rocky planets in disks with varied C/O ratios, as well as h ow the chemica l 
evolution of the galaxy generally can lead to enhanced C/O through time. lLodders! (120041 ) 
suggested that the planetesimals that make up Jupiter's heavy element enrichment were 
carbon- rich, and that Jupiter initially formed at the "tar line" rather than the "ice line." This 
i s one possible expl an ation for the low wa t er abu ndance mea sured by the G alileo Entry Probe 
( jWong et al.l 120041 ) . iKuchner fc Seagerl (120051 ) , similar to IGaidod (120001 ) , were interested 
in giant planets and terrestrial planets that could form in environments where the local 
(or entire disk's) C/0> 1, leading to carbon-dominated (rather than oxygen-dominated) 
silicates. 



More recently IBond et al.l (120101 ) coupled protoplanetary disk abundances derived from 
stellar spectra to a model of disk chemistry, which yields the condensation sequence of solids. 
Their work further c oupled the formation of solids to an N-body model of planet formation 
(jO'Brien et al.l 120061 ). For particular planetary systems, with measured C/O and Mg/Si 
ratios of the host star, they calculated th e equilibrium disk c hemistry and solid composition 
for the initial planetesimal distribution. IBond et al.l (120101 ) furthermore kept track of the 
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contribution of particular planetesimals as they add their mass to growing protoplanets, 
and in the end find the relative contributions of C, O, Mg, Si, etc. to the masses of formed 
planets. 



Within the context of giant planets, iMadhusudhan et al.l ( 1201 lbl ) suggested that day- 
side photometr y of the transiting planet W ASP- 12 b indicates an atmo sphere with C/0> 1.0. 
More recently IMadhusudhan et al.l (j2011al ) and lOberg et al.l (120111 ) have investigated the 
accumulation of gas and icy planetesimals in disks with a range of C/O ratios to understand 
possible pathways to forming "carbon-rich" gas giants. 



1.2. C/O Ratio in Stars 



Composition- dependent planet formation models depend on the stellar ab undances of 



C and O for the initial conditions of disk chemistry . The stellar C/O rati os in iBond et al. 



(hoioh were taken from determinations of C from 



Ecuvillon et al.l (120061 ). Motivate d bv Bond et al.l ( 



201 



01 



cuvillon et al.l (120041 ) and of O from 



larger tabulations of C and O 



abundances were recently made bv lDelgado Mena et al.l (120101) . for 370 FGK stars from the 
HARPS planet-search sample, and IPetigura &: Marcy ( 2011 ). for 457 F and G stars from 
the California Planet Survey sample. These two studies are relatively similar, as they cover 
large samples that include planet-hosting stars and those not found to host planets. 

The two studies do have some differences in the lines of C and O chosen. For the carbon 



abundance, the lDelgado Mena et al.l (120101 ) work used CI lines at 5380.3 and 5052.2 A, with 
only 5380.3 A used for stars with T eS < 5100 K. For oxygen, the forbidden lines of [OI] at 6300 
and 6363 A were used. The derivation of the abundanc es was done w ith a combination of 
the code MOOG, for the generation of synthetic spectra dSnedenl 19731. as updated in 2002), 



the Kurucz ATLAS9 atmosphere grid with overs hooting (iKurucz 



19931 ). and the equivalent 



width s were measured using the ARES program (jSousa et al.ll2007l ). The lPetigura &: Marcy 
( 201ll ) study used a CI line for carbon at 6587 A, and the [OI] line for oxygen at 6300 A. 
The derivation of the abundanc es was performed with the Spectroscopy Made Easy (SME) 
code (jValenti fc Piskunovlll996l ) with Kurucz stellar atmospheres. 



Tabula tions from|Bond et al.l (120101) (who quote Ecuvillon et al. values) . iDelgado Mena et al 



(120 lOi ) , and IPetigura fc Marcyl (120111 ) are shown in Figure [TJ Both of the large studies found 
a somewhat similar sh ape. They found a m aximum at C/O ratios modestly higher that that 
of the Sun f0.55 . fromlAsplund et al.ll2009l) . with a noticeably enhanced peak in the distri- 
bution found bv lDelgado Mena et al.l ( 20101 ). shown most clearly in Figure [T]&. Of particular 
interest to all of these authors, and our Letter, is the tail off to higher C/O ratios > 0.8 
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(dotted line) and even further onto > 1.0 (short dashed line). iDelgado Mena et al.l (12 010) 



find C /0> 0.8 for 24% of their stars, and C/0> 1.0 for 6%. For the iPetigura fc Marcv 
( boilh sample, they find C/0> 0.8 for 29% of their stars, and C/0> 1.0 for 10%. The 
numbers quoted are for the mixed sample of planet-ho sting and non-plan et hosting stars. 
Taken at face value, and the condensation chemistry in iBond et al.l (120101 ). this potentially 
implies carbon planets (formed when C/0> 0.8) in ~25% of planetary systems. 

However, one must take care when estimating the fraction of stars with high C/O ratios, 
given observational error bars. Of interest is the positive tail at high C/O, and a tail such 
as this is expected since the error is approximately constant in the logarithmic abundance 
ratio, [C/O]. This leads to an er ror distribution tha t is lo g- normal in the C/O ratio, and is 
seen in, for instance, Table 6 in IPetigura Marcyl ( 1201 ll ). Their average error at C/0= 1 
is 0.23, which is 1.61<7 above the solar C/O ratio they use, of 0.63. From this 1.61<r, 5.4% 
of the stellar sample is thus expected to be found with C/0> 1, just due to observational 
errors. This would y i eld a "true" fraction of stars with C/0> 1 of ~5%, rather than 10%. 



Delgado Mena et al.l (120101 ) do not provide individual errors on their C/O determinations, 



but assuming similar errors, their entire sample of C/Q > 1 stars (6%) can be explained by 
errors. At C/0= 0.8, in the IPetigura fc Marcyl (120111 ) sample, the average error is 0.16, a 
difference of 0.17 above their solar value, or 1.06a. This yields a n expected fraction of 14.5% , 
which would move the 29% with C/0> 0.8 down to 15%. For bekado Mena et aD fcoioh . 
this expected fraction moves their 24% at C/0> 0.8 down to 10%. The upshot is a much 
reduced tail of stars with high C/O ratios. 



Motivation 



Deriving the abundances of C and O in stellar atmospheres is a dem anding task. Even for 



the Sun, th i s has been an especia lly difficult over the past decade (e.g.. lAUende Prieto et al. 
200 ll . |2002| ; lAsplund et al.l 120091 ) . The most recent state-of-the-art work includes intense 
efforts to find unblended lines, calculated NLTE corrections where applicable, and solar 
atmosp here pressure-temperature profiles that come directly from three-dimensional simu- 
lations (Usplund et al.lbood ; Icaffau et al.lboilh . 



While deriving the abundances of atoms and molecules in low-mass M and K dwarf 
stellar atmospheres is more difficult that than of Sunlike stars, they do have one natural 
asset. Cooler dwarf atmospheres have a natural "flip" in the chemistry at C/0= 1. At 
C/0< 1, most O goes into the CO molecule, but additional O is left over and partitions 
into H2O and, for M stars, TiO/VO gases. TiO/VO absorption bands dominate the optical 
spectra of M stars. However, at C/0> 1, the CO molecules use up nearly all O in cool 
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stars, leaving excess C, which goes into molecules such as C2 and CN. In such a cool star, 
the optical region is dominated by C 2 and CN. We suggest that if C/0> 1 in FGK stars is a 
relatively common phenomenon then cool carbon stars showing C2 and CN bands should be 
relatively common in the sola r neighborhood. This is something that can be investigated, 
and was even briefly noted by iPetigura fc Marcyl (120111 ) . 



True Fraction of Carbon Dwarfs? 



Carbon dwarfs, which appear to have C/0> 1 in their atmosphere, are rare stars. The 
best current understanding of their formation is that they are not made of primordial C/0> 1 
gas, but instead have high C abundances due to accretion of AGB dre dge-up material from a 
companion (e.g., Ide Kool fc Green! Il995l ; ISteinhardt fc Sasselovl 120051 ) . Carbon dwarfs have 
been found by many authors, with mos t being found by large surveys, such as the SDSS 
(IMargon et al.ll2002l ; iDownes et al.ll2004l ). Their numbers are predominantly spread among 
what would otherwise be G, K, and M spectral types. 

For our purposes here, instead of forming these carbon-rich stars by accretion, we can 
take the extremely optimistic view that these are in fact primordial C/0> 1 systems. Then 
the fraction of dwarfs that are carbon dwarfs would be the upper limit on the fraction of stars 
with primordial C/0> 1. Ide Kool &: Green! (119951 ) have previously estimated the frequency 
of carbon dwarfs, based on detections before SDSS. They estimate a space density in the 
disk of ~ 1 x 10 -6 pc— 3, but noted this may be an overestimate. This c an be compared 



to the space density of stars with mass < 0.7 M Q , 6.5 x 10 2 pc— 3, from iBochanski et al. 



(120101 ) . This is a 4-5 order of magnitude difference. 



There are other ways to estimate the carbon dwarf frequency. In particular, ICovey et al. 



(120081 ) have analyzed SDSS spectra and 2MASS photometry of 25,000 sources, in an effort 
to take a census of low mass stars out to J=16.2, over a mass range from 0.1 - 0.7 M . For 
these cool stars, if C/0> 1 did occur, it should be clear, given the abundance of carbon-rich 
molecules (C2, CN) that would occur in such cool atmospheres. Of their sample of 9,649 
low mass stars, they note 24 "exotic contaminants" including 4 carbon stars. This implies 
a C/0> 1 in only 0.04% of cases, although based on color cuts ~l/2 of the carbon dwarfs 
may have been miss ed (K. Covey, personal communication). This work, and the fact that 
Downes et al.l (120041 ) found only ~100 carbon dwarfs from a large SPSS search at 15.6 < 
r < 2 0.8 to specifically find such stars, certainly strengthens the point of Ide Kool &: Green 
(119951 ) that carbon dwarfs are quite rare. The relative frequency is likely 10 -3 — 10 -5 . 



We can think of no reason that Sunlike stars in the solar neighborhood would be uni- 
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formly more enriched in carbon (or oxygen deficient) than the KM stars. One potential way 
out could be if the later-type stars are systematically older than the earlier-type stars, so that 
some amount of galactic chemical evolutio n could have taken place. In this case the younger 
(FG) stars could have higher C/O (e.g. iGaidosI l2000l ). However, the HARPS /California 



samples are not young stars, so constructing a credible explanation through this path seems 
unlikely. 



4. Diagnosing high C/O from Brown Dwarfs? 



The K & M dwarfs from the SDSS give us leverage into the fraction of dwarfs with 
C/0> 1. Are there are any stellar populations that could be used to understand the frac- 
tion of stars with 0.8 <C/0< 1? One possibility are the L and T dwarfs. L-type dwarfs, 
from T e gS of ~ 2400-1400 K, have atmospheres whose infrared opacity is dominated by 
molecular bands of H 2 0, CO, and cloud layers made up of refractory co ndensates like corun- 
dum (AI2O3), enstatite (MgSi0 3 ), forsterite (Mg 2 Si0 4 ), and iron (Fe) (lAckerman fc Marley 
200ll ) . Corundum, enstatite, and forsterite need abundant oxygen to form. If C/0> 0.8, 
or C/0> 1, are common in L-dwarfs, then perhaps we could see a different kind of cloudy 
L-dwarf, with condensates dominated by SiC instead of the Mg-Si-0 forms. Certainly the 
spectra of such objects would be highly abnormal at C/0> 1, as they would strongly fa- 
vor CO at the expense of H 2 0, but perhaps even at C/O ratios closer to 0.8, the cloud 
condensation pathway would be different than normal L dwarfs. 

In the cooler objects, the T-dwarfs, the silicate clouds appear to reside below most of 
the visible atmosphere. Q uite importantly, the se refractory clouds remove 21% of oxygen 
from the atmosphere (e.g. IVisscher et al.l |2010| . and references therei n) . This remo val of 
oxygen from the gas raises the C/O ratio in the visible atmosphere. iLodderj (120101 ) have 
investigated chemistry as a function of C/O ratio at 1100 K and 0.01 bar, which is after 
silicate condensation. Lodders finds a dramatic change in chemistry even at C/0= 0.8. At 
C/0> 0.8, methane becomes more abundant than water (by two orders of magnitude at 
C/0= 0.9) and HCN becomes nearly as abundant as water. Since CO, CH4, H 2 (and 
HCN) all have prominent opacity in the infrared, C/0> 0.8 brown dwarfs could appear 
distinctly different than those with a C/O ratio closer to that of the Sun. 

Since most brown dwarfs are found from color-color diagrams in surveys such as SDSS 
and 2MASS, it is quite possible that high C/O brown dwarfs, should they exist, could have 
previously eluded detection. The cuts from SDSS use i — z color, which could in principle 
be less sensitive to the C/O chemistry since the optical spectra of brown dwarfs are shaped 
by pressure-broadened lines of Na and K ( [Burrows et al.ll2000l ). But 2MASS uses JHK S 
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cuts, which could lead to high C/O dwarfs being missed, since in some cases H2O would not 
be the dominant opacity source, leading to different near-infrared colors. As far as we are 
aware, there are no compelling outliers for high C/O brown dwarfs, but a search for them 
could be valuable. 

A full exploration of C/O ratio in brown dwarfs, and its affects on clouds, is beyond 
the scope of this work. However, there is an additional point that we can motivate. In 
T dwarfs, essentially all infrared opacity is due to water, methane, and carbon dioxide, 
which tie up nearly all C and O. The relatively cloud- free spectra of T dwarfs may allow 
for accurate determinations of C/O in the solar neighborhood. In Figure [2] we show model 
spectra of two 900 K dwarfs, one with C/0= 0.55 (the solar value) and another at C/0= 0.7. 
These models were computed usin g the atmosphere code of M. Marley and J . Fortney and 
collaborators, described elsewhere (IMarley et al.ll2002r iFortney fc Marlevll2007l). The o pacity 



and chemistry databases are taken from a previous tabulation (IFortney et al.l 120051 ) . The 
differences between the two spectra can be prominent, in particular as a ratio. One can 
see differences of 10% in the JHK peaks, and 20-30% in the water bands. With future 
improvements in the accuracy of opacities, in particular a high-temperature database for 
CH 4 , we suggest it may be possible to derive relative C and O abundances, from H 2 0, CH 4 , 
and CO, from infrared T-dwarf spect ra. High spectr al resolution observations of brown 
dwarfs are now being achieved as well (IRice et al.ll2010l ). which could allow further progress 
in abundance analysis. The utility of such efforts would be a better understand i ng of the 
C/O ratio of stars in the solar neighborhood. Recently, some work by lTsuji et al.l (j201lf ) on 
a small number of brown dwarfs with AKARI observations has moved in this direction. 



5. Discussion 

We suggest that recent high-resolution spectroscopic analyses have overestimated the 



Gaidos 


(2000) and 


Bond et al. 


(2010) 



C/0> 0.8 systems will differ in these carbon-rich system, making carbon-rich planets. While 
this logic seems secure, we find it is less likely than some have anticipated. This is because the 
true fraction of carbon-rich parent stars is quite low, 10~ 3 — 10~ 5 . However, our work should 
not be interpreted as claiming that carbon-rich terrestrial or giant planets cannot form. The 
C/O ratio surely varies in the ISM and the region of phase space between 0.8 < C/O < 1 
is less constrained by our work. However, given the overestimation at C/0> 1 it appears 



probable that this region is less-populated than has be en recently suggested (IBond et al. 



20101 ; iDelgado Mena et al.ll2010t iPetigura &: Marcyil201ll ). although not empty. Furthermore 
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we showed in §1.2 that the quoted observational error bars from these groups imply smaller 
measured fractions of high C/O stars, which eliminated 15% from C/0> 0.8 and 5% from 
C/0> 1.0. 

These authors of course put in great effort to understand the sources of their errors, 
and to correct for them. The choice of lines used could be re-examined. The choice of low- 
excitation forbidden [01] lines along with high excitation permitted CI lines could introduce 
systemic errors. A revised (although difficult) study could examine [CI]/ [01] or CI/OI, since 
these lines would behave similarly at a given stellar temperature and in deviations from LTE 
(M. Asplund, personal communication). 



An additional point is that both the large iDelgado Mena et al.l ( 120101 ) and lPetigura fc Marcy 



(120111) works tuned their a bundance retrieval methods to match Sun's C/O ratio, as a stan- 
dard. I Asplund et al.l ( 120091 ). in a recent review, find C/0= 0.55. However, in IDelgado Mena et a 
2010) the particular abundances used for the Sun lead to C/0= 0.66, while in lPetigura fc Marcy 



( 120111 ). it is 0.63. (See the thick tick marks at the top of Figure 1). So there clearly could 
be an offset of ~0.1, towards higher C/O, in these works. 



An issue with the Ni blend of the [OI] line used by iPetigura &: Marcyl (120111 ) may po- 
tentially skew their results towards higher C/O ratios. These authors adopt the abundances 
of O and Ni from 3D abundance analyses, and then use them in their ID analysis. To 
obtain a good fit to the solar spectrum, they adopted log(gf) = —1.98 for the Nil line, 



which is 35% higher th an the laboratory measure ment of -2.11 measured by I Johansson et al. 



( 120031 ). Consequently, IPetigura fc Marcyl 
40% stronger than found in, for example, 
the upward trend in C/O with increasing 
trend is that [Ni/Fe] increases sli ghtly (e.g. 



2011) in their Figur e 3, sh ow a Nil line that is 



Allende Prieto et al.l ( 120011 ). This could explain 



e/H] in their Figure 16, as follows: a known 



metallicity ([Fe/H]) increases (e.g. IDelgado Mena et al.ll2010l ). A Ni blend that is too strong 



Robinson et al.l 120061 ) and [O/Fe] decreases as 



in the solar fit will increase in importance as [Fe/H] increases, leading to an overestimation 
of the Ni blend, an under e stima tion of the O abundance, and a C/O ratio that is too high. 
In the IPetigura fc Marcyl (j201lf ) sample, this could explain why most of the C/0> 1 stars 
occur when [Fe/H]> 0.2. 



We note that IDelgado Mena et al.l (120101 ) and IPetigura fc Marcvl (120 111 ) both relied on 
the stellar a tmosphere models of K urucz. As a comparison, one could alternatively use the 
PHOENIX jHauschildt et ali ll999l) or MARCS JGustafsson et al.ll2008h stellar atmospheres. 
Gustafsson et al.l (120081 ) have compared their MARCS models to published models from 
Kurucz and PHOENIX, and the agreement does indeed vary with T e g and with plane parallel 
vs. spherical symmetry. 



-9 - 



There are still other pathways towards understanding extrasolar abundances. As FG 
stars evolve onto the red giant branch and cool to T c g values of ~ 4000 K, it may be eas- 
ier to constrain their C/O ratios in the same qualitative manner we suggest for M dwarfs. 
The composition of extrasolar planetesimals can be determined by studying externally pol- 
luted white d warf atmos pheres, and some early evidence has been found for carbon-poor 
planetesimals (jJurall2006l ). 



While we have focused exclusively on stellar abundances, it is important to recall that 
condensation of solids within a disk itself can potentially lead to non-stellar C/O ratios in 
nebula gas. For instance, condensation of water in a protoplanetary disk can leave the sur- 
rounding nebula gas relatively carbon-rich (through large abundance o f gaseous CO), which 



1988; 


Lodders 


2010; 


Obere et al. 


2011) 



This may be a pathway to forming giant planets 
with relatively high C/O ratios in their envelopes. A giant planet atmosphere with C/0> 1 
was recently suggest ed for planet WASP-12b, ba sed on fits to 7 photometric points of day- 



side planet emission (IMadhusudhan et al.ll2011bl ). Their finding could be made more robust 



with the inclusion into their model of the absorp tion bands of molecules that are expected 



in the high C/O regime, su c h as H CN and C2H2 (lLoddersll2010l ). which were not considered 



in the IMadhusudhan et al.l (j2011bl ) study. 



Carbon-dominated rocky planets would be extremely interesting objects. Their preva- 
lence around stars of Sunlike abundances, and those stars with enhanced C/O ratios, would 
certainly tell us much about nebular condensation chemistry and the planet formation pro- 
cess. While these planets may be inherently rare, we look forward to additional advances in 
the future. 
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Fig. 1. — Top: Histo gram of the C/O ratios from the papers o f lEcuyillon et al.l (12 004. 
l200fih (as tabulated in lBond et aJboidh , belgado Mena et all told ). and lPetigura fc Marcv 
j201lh . in blue, green, and orange, respectively. C/O ratios of 0.55 (the solar value from 
Asplund et al. 2009), 0.8, and 1.0 are shown as long d ashed, dotted, and s h ort d ashed 
lines, respectively. The adopted solar C/O ratios from iDelgado Mena et al.l (120101 ) and 



Petigura fc Marcyl (120111 ) are shown as thick green and orange ticks at the top. Bottom: 



The data sets a re normalized, yielding frac tions of each sample, in s tead o f number. The 
large samples of belgado Mena et all (120101 ) and IPetigura fc Marcv! (120111 ) find 24-29% of 
FGK stars have C/0> 0.8, and 6-10% have C/0> 1. 
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Fig. 2. — Top: Emitted spectra from two models of a brown dwarf with T e g= 900K and log 
g = 5.0. The thick black curve is for solar abundances, at C/0= 0.55. The thin red curve 
is for a model with C/0= 0.7, with the carbon abundance enhanced. Bottom: The ratio of 
the flux between these two models. 



